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Abstract 

In pattern forming systems such as Rayleigh-Benard convection or direc- 
tional solidification, a large number of linearly stable, patterned steady states 
exist when the basic, simple steady state is unstable. Which of these steady 
states will be realized in a given experiment appears to depend on unobserv- 
able details of the system's initial conditions. We show, however, that weak, 
Gaussian white noise drives such a system toward a preferred wave number 
which depends only on the system parameters and is independent of initial 
conditions. We give a prescription for calculating this wave number, analyti- 
cally near the onset of instability and numerically otherwise. 
47.54.+r, 47.20.Hw, 05.40.-hj, 02.50.Ey 



1 



Typeset using REVT^ 



The classic problem of pattern selection is that of predicting which of a large number of 
available steady states a system will ultimately reach under given experimental conditions. 
In a typical example, that of directional solidification, there is a simple steady state of the 
system in which the solidification front is planar and advances into the melt at a constant 
speed. By varying a control parameter, one reaches a regime in which this state is linearly 
unstable against all perturbations with wave numbers in a given interval. In this regime 
the front settles into a spatially periodic, cellular shape. Ample evidence exists that there 
is one such patterned steady state with each wave number q in the interval of instability. 
Moreover, there is a finite subrange of wave numbers for which the cellular steady states are 
themselves linearly stable. The pattern selection question is then this: Into which of these 
patterned states will the system restabilize in a given experiment? 

The answer to this question appears to be that the final wave number depends not only 
on the system parameters, but also on the details of the initial conditions from which the 
system evolves. Since these details cannot be observed in practice, the final wave number 
is not reproducible. We will argue in this Letter, however, that among the possible steady 
state wave numbers there is one which is preferred, in the sense that subjecting the system 
to weak, Gaussian white noise drives it toward that wave number, and in the long-time limit 
that wave number is overwhelmingly more probable than any other. We will show how this 
preferred wave number can be calculated for one-dimensional systems. 

Noise effects on pattern-forming systems have been studied by many authors, although 
few have considered the role of noise in readjusting the wave number of an established, 
periodic pattern. Deterministic evolution from random initial conditions has been studied 
in amplitude equations and noise effects on the initial stages of pattern formation 

has been investigated in Benard convection |Q, dendritic growth and cellular ^ and 
dendritic |^ arrays in directional solidification. Several numerical studies of the Swift- 
Hohenberg equation (which is relaxational) have all shown p|-p!3[] that noise selects the wave 
number which minimizes the underlying free energy, as expected. In a series of papers 
which are particularly relevant to our work, Kerszberg carried out numerical simulations of 
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directional soldification with and without Gaussian white noise [0,0. This is a system 



which has no underlying free energy. With noise added, the system restabilized into a 
cellular state with a unique, reproducible wave number independent of initial conditions, 
while without noise the final wave number depended on initial conditions. 

The relaxational, or "gradient", case is quite helpful for motivating our calculations. 
Suppose the state of the system is specified by giving a set of amplitudes Xk of Fourier 
modes of wave number A;, chosen so that the simple steady state of the system is Xk = 0, 
and that the control parameter has been set so that this state is linearly unstable. (We will 
also take the Xk to be real, so that the pattern is left-right symmetric.) The amplitudes 
evolve according to 

dxk d^{x) 
dt dxk 

where $(2;) is the free energy. This evolution always makes $ decrease with time. A 
patterned steady state with wave number q will have Xk nonzero only when k is an integer 
multiple of g; let the value of the free energy of this state be (^ss{<l)- If we add noise of strength 
e to (|I]), it will produce occasional large fluctuations which take the system far enough out 
of the local free energy minimum at the state of wave number q that it then relaxes to a 
different local minimum, with wave number q' . The relative probability of a transition from 
q to q' versus a transition from q' back to q is proportional to exp[($s5(g) — ^ss{q'))/e\. That 
is, fluctuations are more likely to take the system from a state of higher free energy to one 
of lower free energy than vice versa. In pattern forming systems, the steady states near 
the edges of the band of stable wave numbers, which are almost unstable, have higher free 
energies than states in the interior of the band, so noise drives the system away from those 
states and toward the absolute minimum of the free energy. 

We argue here that the same thing happens in non-gradient systems, i.e., ones in which 
the time evolution is given by 
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but the Fk cannot be put into the form (|l]). To see this, we first add noise to the dynamics 
to obtain the Langevin equation 



dt 



(3) 



where the are independent Gaussian random variables with mean zero and unit variance, 



(4) 



The Xk are then random variables, and a standard argument [|l^] leads to the Fokker-Planck 
equation, which describes the time evolution of their probability distribution. In the long- 
time limit, this converges to a steady-state distribution Vss{x), which satisfies 



d 
dxk 



e dV 

Fk{x)Vss{x] 



0. 



(5) 



2 dxk 

Here and below, repeated indices are to be summed. Since the noise terms in (0) are 
independent of x, there is no difference here between the Ito and Stratonovich interpretations. 

In the weak-noise limit e — 0, we can solve using a WKB method Since 
probability distributions must be positive, we may write 



Vss{x) = exp[-S'(a;)/e] 



(6) 



so that S plays a role similar to that of the free energy in the gradient case; in fact the 
solution of is S' = 2$ in that case. From (|]) we then find that to leading order in e, 5" 
satisfies 



dS 
dxk 



1 dS 



2 dxi 



+ Fkix) 



0. 



(7) 



From this we see that any steady state of the deterministic dynamics (^, i.e., a point where 
all Fk vanish, is a stationary point of 5*. We will see below that a linearly stable steady state 
of is in fact a local minimum of 5*. 

We can also see that S plays the role of $ in a different sense. The deterministic evolution 
(S) never makes S increase, since 



— - —F (x) - — — <0 (S) 
dt dxk ^ 2 dxk dxk ~ ' 

Thus even a non-gradient system can be considered relaxational, since there is a function 
S{x) which is non-increasing. In gradient systems, however, that function is just the free 
energy, which is generally easier to compute. 

Eqn. can be solved by the method of characteristics, or equivalently by recognizing 
it as the Hamilton- Jacobi equation, with energy 0, for the Hamiltonian 

H{p, x) = ~PkPk - PkFk{x). (9) 

The characteristic curves are given by 

ik = -Pk - Fk{x), Pk = Pk'^^, (10) 

dxk 

where overdots represent derivatives with respect to some parameter r. Along a character- 
istic, S is given by 

S = Pkik = -Pk{Pk + Fk). (11) 

Since we have TC = 0, this becomes S = —{pkPk)/2, so 5* always decreases along a charac- 
teristic. 

To find the preferred wave number, we need only calculate S{x) at each stable steady state 
of (^. The wave number of the state with the lowest value of S is the preferred wave number: 
by (^), if the wave numbers are discrete then all others are exponentially less probable. 
However, since S always decreases along a characteristic, there are no characteristics which 
run from one local minimum of S to another. Thus in order to compare the values of S at 
its various local minima, we need characteristics running from one common point to each 
minimum. Fortunately, we can find these - the common point is the unstable steady state 
Xfc = 0, the simple steady state of the system. That there should be a characteristic linking 



this to each stable steady state comes from a counting argument. Fixed points of (|TOD occur 
where = and Fk{x) = 0; linearizing about any fixed point gives 



Sx = -p- MSx, p = M'^p, (12) 



where M is the matrix which gives the hnear stabihty of the deterministic dynamics 



(13) 



dxk' 

If there are amphtudes x^, then the phase space {x,p) is 2A^-dimensional; from ( [I^ ) we 
see that there are N eigenvectors at each fixed point with eigenvalues which are the same 
as the stabihty eigenvalues of the corresponding deterministic steady state, and another 
eigenvectors with the negatives of those eigenvalues - and all p components equal to 0. From 
(p!0|), we see that the p = subspace of phase space is invariant, and ([lID shows that following 
a trajectory in this subspace would leave S unchanged. Thus we want characteristics not 
to be in this subspace. However, the only eigenvectors coming out from a linearly stable 
steady state (i.e., having positive eigenvalues) are those which lie in the p = subspace. 
Thus we cannot have a characteristic leaving a stable steady state. All nearby characteristics 
must then be directed toward that state. Since S always decreases along characteristics, a 
linearly stable steady state of the deterministic dynamics must then be a local minimum of 
S. However, the steady state at x = is not stable; it has one linearly unstable direction 
for each fundamental wave number q against which it is unstable. Correspondingly, for 
each such q there is a steady state which bifurcated from the x = state as the control 
parameter was increased to its present value. It is then natural to expect that there will 
be one characteristic running from {x,p) = (0,0) to each linearly stable steady state (and 
p = 0). 

To do the calculation numerically, it is useful to have the expression for the Lagrangian 
corresponding to the Hamiltonian (^): 

C{x,x) = -^[xk + Fk{x)][xk + Fk{x)] (14) 

From (10), we see that on the characteristics this is equal to —{pkPk)/'2, which in turn is equal 
to S. Thus we can set up the calculation of S at the steady state as a minimization problem, 
using Hamilton's principle: we wish to find the trajectory which approaches {x,p) = (0, 0) as 
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r — oo, approaches the steady state x (and p = 0) as r — > oo, and minimizes / Cdr. The 
resulting minimum value is the value of 5* for that particular steady state. Minimizing this 
result in turn over the possible fundamental wave numbers q gives the preferred wave number. 
This procedure is made much easier by the fact that the structure of the deterministic 
equations is usually such that the multiples of any fundamental wave number q form an 
invariant subspace. This property is inherited by the equations (p!oD for the characteristics, 
and so only a relatively small number of modes must be kept for each q. 

When the control parameter is only slightly beyond the onset of instability of the simple 
steady state, it is possible to calculate the preferred wave number analytically. The simplest 
nontrivial example of this procedure is for a system in which a symmetry keeps the k = 
amplitude fixed. Expanding to third order in the interface displacement generally leads to 
equations of the form (H) with 

F2q = -\cr2q\x2q + CXX^, 
'3q = -\03q\X3q f^^g 



03q I X'iq + (3x1 + ^^q^'^q ■ (15) 



All coefficients depend on the fundamental wave number g, and near the onset of instability 
of the state a; = the linear growth rate aq is small and positive. The p^q equation coming 
from (|10D immediately gives p^q oc exp(— |cr3g|r), and since p^q must go to zero for r ^ — oo 
this forces p^q = 0. This makes x^q irrelevant, since it does not appear in Fq or F2q, and in 
5* it is multiplied by p^q. The patterned steady state is given by 



<= Wq\a2q\/{\\a2q\+afi)Y/' + 0{a'J' 
4q= o-g"/(A|cr2g| + afi) + 0{aq 



xZ= aqa/{\\(y2q\ + a/i) + ©(aj), (16) 



with the corrections coming from higher order terms which were neglected in writing Eqn. 
(p!5|). With the rescalings 



Xq = ^y^a;, Pq = af^p, X2q = (TqV, P2q = Cr'^T, (17) 



we find that x and p are explicitly of order cTq, while to leading order we always have 

y = ax'^/\a2q\, r = -fip/\a2q\. (18) 

Substituting all this into the expression (|]) for the Hamiltonian and setting it equal to zero 
then gives 

p= -2x + 2{X\a2q\+afj,)xy\a2g\. (19) 

It is then simple to integrate (0) to leading order to get the value of S at the steady state: 

S{x'') = -aJ|a2,|/2(A|a2,| + a/x). (20) 

To find the preferred wave number q to leading order, we must minimize (^) over all q in the 
unstable band. Higher-order corrections can also be calculated, although to do this properly 
we need to retain higher-order corrections in (pTS]), which in turn requires us to keep more 
than three Fourier modes. Note that in the special case fi = —2a, ( [T5| ) is a gradient system 
(to the relevant order in ag). However, this leads to no simplification of the calculation of 
S, although the result then agrees with the free energy. 



We have carried out the analogous calculations for the Greenside-Cross equation for 
Benard convection, which is a non-gradient equation for a spatial order parameter ip, 

d^/dt = [7 - (V^ + + 3\Vij\^V^t/j. (21) 

The steady state ip = is linearly unstable for 7 > 0. Keeping modes up to the seventh 
harmonic, we find that the preferred wave number for small 7 is given by 

, = i_2_12l2!_£5l2!_oM. (22) 

^ 4 1024 16384 ^ ' ^ ^ ^ 

In one spatial dimension, the Greenside-Cross equation is a gradient system, and this result 



agrees with the expansion obtained by minimizing the free energy. Cross and Meiron [jl3 
in numerical simulations of deterministic, two-dimensional evolution with random initial 
conditions, find that the system reaches a wave number q ^ 0.78 for long times for 7 = 1/2. 



Setting 7 = 1/2 in ( ]22|) gives q = 0.843, although the series has clearly not converged (the 
first three terms give 0.850). 

In this Letter we have argued that pattern forming systems have a naturally preferred 
wave number, namely that which the system would approach if it were subjected to weak, 
additive, Markovian, Gaussian white noise. This is not to say that the noise actually ex- 
perienced by such a system has any of these five properties; but Gaussian white noise is 
appropriate because it does not bias the system's preference for a wave number, as colored 
noise, for example, would do. The argument does suggest, however, that the role of noise 
is not limited to providing the initial fluctuation which takes the system out of its simple 
steady state and starts it evolving toward one of its possible restabilized states, as is usually 
taken (implicitly) to be the case. Rather it has a continuing role in readjusting the wave 
number of the patterned state, generally through occasional large fluctuations which either 
create or destroy a cell. 

An obvious question which needs to be addressed in the future is that of the rate at 
which the steady-state probability distribution Vss is established. It should be possible to 
use the same classical-mechanics techniques to study this question as we have used above, 
since substituting the time-dependent generalization of ansatz (^) into the Fokker-Planck 
equation leads directly to the time-dependent Hamilton- Jacobi equation for the Hamiltonian 
(^). Further work on this point is under way. 

Other important issues arise when we model the noise to which a system is actually 
subjected, which may not be additive, white, Gaussian, Markovian, or weak. For instance, 
the relevant fluctuations might be in the value of the control parameter. If any of the first 
four properties are lacking, then the appropriate Fokker-Planck equation will not have the 
simple form (|^). Even if the relevant noise source is thermal, evolution equations of the form 
(0) often arise only after considerable manipulation of some more complex model which is 
written down from first principles. The e in (|^) may then be replaced by something which 
depends on q and even x. In all such cases the subsequent calculations need to be modified 
appropriately. 
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When the noise strength is finite, several important effects arise. One is that we may 
need the next higher order correction to the leading-order S{x) which we have calculated. 
A second, related point is that the relevant probabilities to compare are not just the heights 
of the peaks in Vss, but the areas under the peaks - including not just the exact steady 
states, but also perturbations of those states. These effects have been seen by Kerszberg []T^, 
who found that even in a gradient system the observed wave number was not equal to the 
wave number which minimized the free energy when noise was included in his calculations. 
Finally, in an infinite system the possible band of wave numbers is continuous, so that it 
is not true that one wave number has a probability which is exponentially larger than all 
others. Rather the probability exp[— S'(g)/e] is appreciable for all wave numbers within a 
range of order e^^^ around the "preferred" wave number. Thus if the noise variance is larger 
than something of the order of the inverse square of the size of the system, then there is still 
a range of possible wave numbers that one might observe even in the long time limit, albeit 
a very narrow range for finite but weak noise. 

Finally, it would of course be far more satisfying to have a direct physical interpretation 
of the criterion for the preferred wave number, rather than just a prescription for calculating 
it. 

This work was supported by the National Aeronautics and Space Administration, 
through grant number NAG3-1603 and also through the JOVE program at North Dakota 
State University. 
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